ADA103835 


Source  Flow  Effect  on  Lineshape 


H.  MIRELS 

Aerophysics  Laboratory 
Laboratory  Operations 
The  Aerospace  Corporation 
El  Segundo,  Calif.  90245 


10  July  1981 


* 


APPROVED  FOR  PUBLIC  RELEASE; 
DISTRIBUTION  UNLIMITED 


Sponsored  by 

DEFENSE  ADVANCED  RESEARCH  PROJECTS  AGENCY  (DOD) 

DARPA  Order  No.  3646 

Monitored  by  SD  under  Contract  No.  F04701-80-C-0081 

SPACE  DIVISION 
AIR  FORCE  SYSTEMS  COMMAND 
Los  Angeles  Air  Force  Station 
P.O.  Box  92960,  Worldway  Postal  Center 
Los  Angeles,  Calif.  90009 

D 


DTIC 


ELECTE 
SEP  8  1981 


D 


THE  VIEWS  AND  CONCLUSIONS  CONTAINED  IN  THIS  DOCUMENT  ARE  THOSE 
OF  THE  AUTHORS  AND  SHOULD  NOT  BE  INTERPRETED  AS  NECESSARILY 
REPRESENTING  THE  OFFICIAL  POLICIES,  EITHER  EXPRESSED  OR  IMPLIED,  OF 
THE  DEFENSE  ADVANCED  RESEARCH  PROJECTS  AGENCY  OR  THE  U.S 
GOVERNMENT. 


3 1  9 


This  report  was  submitted  by  The  Aerospace  Corporation,  El  Segundo,  CA 
9024S,  under  Contract  No.  P04701-79"e*00W  with  the  Space  Division,  Deputy  for 
Technology,  P.0.  Box  92960,  Worldway  Postal  Center,  Los  Angeles,  CA  90009.  It 
was  reviewed  and  approved  for  The  Aerospace  Corporation  by  W.  R.  Warren,  Jr., 
Director,  Aerophysics  Laboratory,  Lieutenant  J.  C.  Garcia,  SD/YLVS,  was  the 
project  officer  for  Technology.  Dr.  H.  Allan  Pike  is  the  Program  Director  for 
the  DARPA  Washington  Office.  This  research  was  supported  by  the  Defense 
Advanced  Research  Projects  Agency  of  the  Department  of  Defense. 

This  report  has  been  reviewed  by  the  Public  Affairs  Office  (PAS)  and  is 
releasable  to  the  National  Technical  Information  Service  (NTIS).  At  NTIS,  it 
will  be  available  to  the  general  public,  including  foreign  nations. 


This  technical  report  has  been  reviewed  and  is  approved  for  publica¬ 
tion.  Publication  of  this  report  does  not  constitute  Air  Force  approval  of 
the  report's  findings  or  conclusions.  It  is  published  only  for  the  exchange 
and  stimulation  of  ideas. 


Jdme 8  C^e^frc ia/  Lt,  OSAF 
Project  Officer 


Florian  P.  Meinhardt,  Lt  Col,  USAF 
Director  of  Advanced  Space  Development 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  (-When  Data  Enfrtd) 


naroi  W fj 4 . 1  / 

■  Entered)  [ 


REPORT  DOCUMENTATION  PAGE 


*•  TITLE  (and  Subtitle) 


2.  GOVT  ACCESSION  NO 

2  %D'£tf3  $35 1 


['SOURCE  FLOW  EFFECT  ON  LINESHAPE 


7  AuTHORf.J 

Harold  Mirels 


/ 


READ  INSTRUCTIONS 
DEFORE  COMPLETING  FORM 


3.  RECIPIENT'S  CATALOG  NUMBER 


5  TYPE  OF  REPORT  A  PERIOO  COVERED 


[/*•  performing  owe:  report  number 
k  TR-008 1(6764)-! 


I*.  CONTRACT  OR  GRANT  NUMBERfa) 

^^F04701-80-C-008L./ 


|9  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

J The  Aerospace  Corporation 
El  Segundo,  Calif.  90245 


II.  CONTROLLING  OFFICE  NAME  AND  ADDRESS 

Defense  Advanced  Research  Projects  Agency 
1400  Wilson  Blvd. 

Arlington,  VA  22209 


// 


U  MONITORING  AGENCY  NAME  A  ADDRESSf/f  dllleron-  trom  Controlling  Ollict) 

Space  Division 

Air  Force  Systems  Command 

Los  Angeles,  Calif.  90009 


IS.  DISTRIBUTION  STATEMENT  (of  thtm  Roport) 

Approved  for  public  release;  distribution  unlimited. 


program  element,  project,  task 

AREA  A  WORK  UNIT  NUMBERS 


»2.  REPORT  DATE 

10  July, 1981 


13.  NUMBER  OF  pages 

15 


* 


15.  SECURITY  CLAS5T  (ofth! a  report) 

Unclassif ied 


ISa.  DECLASSIFICATION  DOWNGRADING 
SCHEDULE 


17.  DISTRIBUTION  STATEMENT  (ol  the  abatract  entered  In  Block  20,  II  different  from  Report) 


18.  supplementary  NOTES 


|  19.  KEY  WORDS  ( Continue  on  reverae  aide  If  neceaaary  and  Identity  by  block  number) 

CW  Chemical  Laser 
Inhomogeneous  Broadening 
Source  Flow 
Spectral  Lineshape 


20.  ABSTRACT  (Continue  on  reverae  aide  If  neceaaary  and  Identify  by  block  number) 

The  effect  of  mean  fluid  motion  on  the  lineshape  of  a  Doppler -broadened  medium 
is  investigated  for  the  case  of  radiation  that  is  perpendicular  (transverse)  to 
the  axis  of  a  two-dimensional  source  flow  with  semi  angle  8  .  The  case  <  <  1 
is  considered.  Transverse  flow  effects  are  significant  when  (V  /a)2  >0(f),  where 
Ve  and  a  characterize  the  transverse  mean  motion  and  thermal  motion,  respectively 
For  typical  cw  chemical  lasers,  the  latter  condition  corresponds  to  flows  in  the 
range  0£  >  0  (1/25). 


00  F0RM  1473 

(FACSlMILd 


UNCLASSIFIED 


/ 


,  SECURITY  CLASSIFICATION  OF  THIS  PAGE  (Whmn  Dim  Bntorod) 


CONTENTS 


I.  INTRODUCTION . 5 

II.  THEORY . 7 

III.  APPLICATION . 17 


A'crssion  For 
|  F'tt'  C.FA&I 
DTIC  TB 
Unanncireed 
Just if iceti"n 


Vtxr 

— - - ■ — - 

Distribution/ 

Availability  Codes 

1 - 

Avail  and/or 

Dist 

fl 

Spec 

Lai 

DTIC 

SELECTE 
SEP  8  1981 

D 


i 


FIGURES 


1.  Source  Flow . 6 

2.  Particle  Space  in  V-v  and  W-w  Coordinates.......... . 10 

3.  Llneshape . 12 

4.  Anomalous  Index  of  Refraction . 15 

5.  Effect  of  Threshold  Gain  and  Source  Flow  on  Output  Power  from 

Multiple  Longitudinal  Mode  CW  Chemical  Laser  Employing  a  Fabry- 
Perot  Resonator . . . . . . . 18 


3 


I .  INTRODUCTION 


The  effect  of  particle  motion  on  spectral  lineshape  (i.e.t  Doppler 

broadening)  is  usually  evaluated  by  assuming  that  the  particles  have  a  random 

e 

thermal  motion.^  In  high-speed  gas-flow  lasers,  of  which  chemical  lasers  are 

'jy 

an  example,  the  working  fluid  may  have  significant  mean  velocities  as  well  as 
random  velocities  in  the  optical  path  direction.  In  these  cases,  it  is  neces¬ 
sary  to  take  the  mean  motion  into  account  when  evaluating  the  spectral  line- 
shape.  The  effect  of  a  source  flow  on  Doppler-broadened  lineshape  is  eval¬ 
uated  herein  for  the  case  where  radiation  is  perpendicular  to  the  source  flow 


axis 


I  — 

,ClTgr~i7. 


The  application  to  a  cw  chemical  laser  is  then  noted. 


f 


*A.  Yariv,  Introduction  to  Optical  Electronics,  (Holt,  Rinehart  and  Winston, 
New  York,  19^1)  pp.  76-79. 

^D.  J.  Spencer,  H.  Mirels,  and  D.  A.  Durran,  J.  Appl.  Phys.  43,  1151  (1972). 
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II.  THEORY 


Consider  a  two-dimensional  source  flow  with  half  angle  9g  (Fig.  1).  A 
radiation  field  I(v)  is  assumed  to  propagate  in  a  direction  that  is  transverse 
to  the  flow  axis.  The  divergence  of  the  flow  results  in  mean  motion  in  the 
transverse  (optical  path)  direction.  Local  mean  velocity  and  random  velocity 
in  the  transverse  direction  are  denoted  by  V  and  v,  respectively.  The  resul¬ 
tant  transverse  velocity  field  is  W  =  V  +  v.  The  radiation  frequency  v,  which 
is  resonant  with  particles  of  velocity  W,  is  found  from  the  Doppler  relation^ 

W  =  (v  -  vq)  c/vo  (1) 

where  vq  is  the  resonant  frequency  for  stationary  particles  and  c  is  the  speed 
of  light.  The  spectral  lineshape  is  found  by  determining  the  distribution 
function  for  particles  in  the  range  W  to  W  +  dW,  which  is  found  in  the  fol¬ 
lowing  paragraphs. 

2 

Consider  the  case  9g  «  1.  The  flow  density  and  axial  flow  velocity  at 
each  streamwise  station  x  are  independent  of  9,  whereas  the  transverse 
velocity  can  be  expressed 


(2) 


where  Vfi  is  the  value  of  V  corresponding  to  0  ■  9e-  The  fraction  of  particles 
in  the  velocity  range  V  to  V  +  dV  is 


F(V)  -  1/(2  Vfi)  IVl  <  Ve 

-  o  I vl  >  ve  (3) 

Under  equilibrium  conditions,  the  fraction  of  particles  in  the  velocity  range 
v  to  v  +  dv  has  a  Maxwellian  distribution* 

i  2 ,2 

f(v>  -  -i-  «-V  /a  (4) 

ir  a 

where  a  =  (2  kT/m)*^  is  the  most  probable  random  particle  speed.  The 
quantities  F(V)  and  f(v)  are  normalized  so  that 


r 

r 


f(v)dv  =  r°  F(V)  dv  =  1 


f ( v)  F(V)dv  dV 


1 


(5a) 


(5b) 


The  integrand  in  Eq.  (5b)  denotes  the  fraction  of  particles  in  the  combined 
range  V  to  V  +  dV  and  v  to  v  +  dv.  We  now  replace  the  independent  variables 
V,  v  by  the  pair  W,  w  defined  by 


W  -  v  +  V 


(6a) 


8 


w 


v  -  V 


(bb) 


It  is  seen  from  Fig.  2  that  W,  w  are  orthogonal  variables.  The  particles  are 

confined  to  the  region  -V  <  V  <  V  in  V,  v  space  and  to  the  region  W  -  2  V 

<  «  <  H  +  2  Vg  in  M,  w  space.  The  particles  in  the  velocity  range  V  to  V  +  dV 

and  v  to  v  +  dv  are  related  to  the  particles  in  the  range  W  to  W  +  dW  and  w  to 

w  +  dw  by 


f(v)  F(V)  dvdV  =  f(—  )  F(~) 


3(v,V) 

3(w,W) 


dw  dW 


(7) 


O 

where  3(v, V)/3(w,W)  is  the  Jacobian  given  by 


3(v,V)  =  iVw  Vw,  „1 
3(w,W)  'vw  Vw'  2 


Equation  (5b)  becomes 


/“  /»W+2 

dW  J 

—  W-2V 


W+2Ve  -(vrW)2/(4a2) 

dw  -  - 

4ir  '  aV 


) 


£ du  [ 


w+v  w-v 

erf( - e)  -  erf( - e) 

K  a  J  K  a  1 1 

4V 


(8) 


(9a) 


(9b) 


The  integrand  in  Eq.  (9b)  can  be  interpreted  as  the  fraction  of  particles  in 
the  range  W  to  W  +  dW.  Introduce  the  notation  Xg  =  Vg/a,  and 

•^D.  H.  Menzel,  Fundamental  Formulas  of  Physics,  Vol.  I,  (Dover  Publications, 
New  York,  1960) ,  p.  26^ 
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(10) 


..  V-V  ,  /0  v-v 

X  =  -  =  - 2  -c-  =  2(ln  2) 1/2  ° 

a  v  a  v  ' 


Av„ 


where  Av^  is  the  Doppler  width  (FWHM)  corresponding  to  Xg  =  0.  The  fraction 
of  particles  in  the  interval  X  to  X  +  dX  is  denoted  n(X,  Xg)  and,  from  Eq. 
(9b),  equals 


n(X,  X  )  =  [erf  (X  +  Xe)  -  erf(X  -  Xe)]/(4X  )  (11) 


Limiting  forms  of  n(X,Xg)  are 


(12a) 


-  e'X  [1  +  0(X2)] 


(12b) 


,l/2 

2X 


erf  (Xj[l  +  O(X^) ] 


(12c) 


2 

Thus,  transverse  mean  flow  effects  are  negligible  for  (Xg)  <<  1  and  must  be 
considered  for  (Xfi)2  >  0(1).  This  result  is  physically  realistic,  since  Xe  is 
the  ratio  of  characteristic  mean  to  characteristic  random  motion  in  the  trans¬ 
verse  direction.  The  variation  of  it 1/2  n(X,Xe)  with  X  is  plotted  in  Fig.  3  for 
various  values  of  Xg.  Figure  3  can  be  used  directly  to  estimate  the  effect  of 


11 


transverse  mean  motion  on  zero  power  gain  of  a  Doppler-broadened  medium.  If 


we  let  gQ(X,  Xe)  denote  the  zero  power  gain  corresponding  to  X,  Xe,  it  follows 

that  in  the  Doppler  limit4*5  (i.e.,  in  the  limit  of  Av /Av  <  <  1,  where  Av. 

n  D  n 

is  the  homogeneous  half-width) 


g0(x,xe) 

go'(0'J0 


X1'2 


n(x,Xe) 


(13) 


where  go(0,0)  is  the  line  center  value  in  the  absence  of  transverse  motion. 
Thus,  the  ordinate  in  Fig.  3  is  a  direct  measure  of  the  effect  of  transverse 
mean  motion  on  zero  power  gain.  Line  center  gain  is  decreased,  and  the 
lineshape  is  broadened  as  X_  is  increased.  The  corresponding  value  of  the 
anomalous  index  of  refraction  n(X,Xe)  can  be  expressed,4,5  for  X  =  c/vq. 


2. 

X  ■  g(0,0) 


dX 

rf  [erf(Xo+Xe)  -erf(Xo-xJ] 

O 


■T72  r  “M1  +  °(x2l) 


(14a) 


(lAb) 


-  — /2D(X)[1  +  0(x*)]  (14c) 


V  Mi re Is,  AIAA  J.  (5),  478  (1979) 

5M.  Sargent,  M.  0.  Scully  II,  and  W.  E.  Lamb,  Jr.,  Laser  Physics,  (Addison- 
Wesley,  Massachusetts,  1974),  pp.  144-155. 
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where  D(  )  is  the  Dawson  integral.  Equation  (14a)  is  plotted  in  Fig.  4.  The 
maximum  value  of  the  index  decreases  as  Xg  increases  and  occurs  at  values  of  X 
and  Xe  related  by  D(X  +  XQ)  -  D(X  -  Xe),  which  for  (X  +  Xe)“2  «  1  becomes  X2 

=  x2e  +  0.5  [i  +  o(x  +  xer2]. 
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III.  APPLICATION 


It  is  convenient  to  express  the  parameter  Xg  in  terms  of  local  flow  Mach 
number.  The  most  probable  speed  a  is  related  to  the  local  speed  of  sound  a8 
-  (yRT/m)1^  by  the  relation  ag/a  ■  (y/2)*^2  where  y  is  the  ratio  of  specific 
heats,  which  equals  7/5  and  5/3  for  diatomic  and  monatomic  gases,  respectively. 
Let  Mg  5  Ve/ag  denote  the  characteristic  transverse  flow  Mach  number,  and  it 
follows  that  Xe  is  related  to  Mg  by 


Xe  -  (2/y)1/2  Me  (15) 

Thus,  Xe  is  nearly  equal  to  the  characteristic  transverse  flow  Mach  number. 

The  present  results  can  be  applied  to  cw  chemical  lasers  in  which  trans- 
verse  flow  expansion  is  permitted.  Typical  axial  flow  Mach  numbers^  are  in 
the  range  4  to  6  and  can  be  characterized  as  being  of  order  5.  The  cor¬ 
responding  value  of  Xe  is  then  of  order  Xe  ■  0(50^).  Thus,  transverse  flow 

effects  are  negligible  for  02  «  0(1/25)  and  must  be  considered  for  0^  >  0(1/25). 

e  e 

The  effect  of  source  flow  on  cw  chemical  laser  output  power  is  being  inves¬ 
tigated  for  the  case  of  a  two-level  model, ^  laminar  mixing,  a  Fabry  Perot 
resonator,  and  multiple  longitudinal  modes  spacing  [i.e.,  Av^/Av^  <0  (1), 
where  Av^  ■  c/2L  is  the  longitudinal  mode  spacing].  Preliminary  results  are 
indicated  in  Fig.  5.  The  ordinate  P/Psat  is  the  ratio  of  output  power  to 
saturated  output  power,  and  the  abscissa  is  the  ratio  of  cavity  threshold 
gain,  gc,  to  the  maximum  zero  power  gain,  g^p*  for  the  case  Xg  ■  0.  The 
decrement  in  output  power  caused  by  source  flow  is  seen  to  become  more  severe 
as  threshold  gain  is  increased. 
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